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The optimal design for infrared second harmonic generation (SHG) is determined for a GaAs- 
based quantum device using a recently developed genetic approach. Both compositional parameters 
and electric field are simultaneously optimized, and the quantum limit for SHG, set by the trade- 
off between large dipole moments (favouring electron derealization) and large overlaps (favouring 
electron localization), is determined. Optimal devices are generally obtained with an asymmetric 
double quantum well shape with narrow barriers and a graded region sideways to the largest well. 
An electric field is not found to lead to improved SHG if compositional parameters are optimized. 
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INTRODUCTION 

In semiconductor quantum structures intersubband 
transitions lead to very large second-order optical 
susceptibilitiesa; in particular, asymmetric quantum con- 
finement may induce large second harmonic genera- 
tion (SHG) in material systems, like Si/SiGe, where 
SHG is prohibited by .symmetry in the centro-symmetric 
bulk host materialso; in non-centrosymmctric sys- 
tems, like GaAs/AlGaAs, quantum confinement-induced 
SHG may exceed by orders of magnitude the bulk 
contributiontrO. Intersubband SHG has been intensely 
investigated in these and other classes of materials, like 
AlInAs/GalnAsQla or GaN/AlGaNH; band-jto-band non- 
linear susceptibility has also been studiedl!3il3 

The intersubband SHG coefficient is given bya 
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where Ei m = Ei — E m is the transition energy between 
subbands I and to, p\ m = (l\z\m) is the dipole matrix 
element of the transition, hu> is the photon energy, pi 
is the surface charge density of subband I, and Ti m ac- 
counts for inhomogenous broadening of the / — ► m tran- 
sition. Most studies of SHG focus on the double-resonant 
regime, where TUo is in resonance with two intersubband 
gaps; in this regime two denominators in (|l|) are simul- 
taneously resonant, and the device can be scketched as 
a three level system with the frequency of the incoming 
electromagnetic field matching approximately the 1 — > 2 
and 2 — > 3 transitions. Hence, at the resonant frequency 
3 pi/(e r 2 )^(w r ), where 
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Quantum confinement-induced SHG (i.e., in addition 
to the possible bulk contribution) can only be observed 
if the confinement potential does not have the reflection 
symmetry, so that none of the dipole matrix elements 
in (|^) is zero; thisjpan_be achieved either by an asym- 
metuc-campositionEI & i3 or by application of an electric 
fieldO'Eirila F or bothQ. However, the use of an electric 
field, in addition to a built-in asymmetry, does not neces- 
sarily improve SHG provided that resonance is achieved 
by proper design of the heterostructure, so that in gen- 
eral the usefulness of an applied electric field from this 
point of view has still to be assessed. 

Double-resonant SHG can be maximized by appropri- 
ately engineering the structure, in order to make the 
product of the three dipole matrix elements appearing 
in (||) as large as possible; to this end, the best compro- 
mise between large intersubband dipole moments (which 
favours electron derealization) and large overlaps (which 
favours electron localization) must be found. The opti- 
mal design of the confinement potential, however, is a 
very hard optimization problem, since, in general, many 
compositional and geometrical parameters, and possibly 
the external field, can be varied, while resonance with a 
given radiation should be preserved in the process; fur- 
thermore, the optimal structure depends on the radiation 
frequency for a given class of materials. This optimiza- 
tion task has been undertaken so far only for specific 
confinement profiles; apart from the required asymme- 
try, the selected shapes were simple enough that opti- 
mization could be perfapiijed with respect to a single 
parametert! or the fieldEHl3'll3. Analytical xuitimization 
methods for SHG have also been proposedE0E3, provid- 
ing a qualitative estimate of the optimal SHG which can 
be obtained in a heterostructure; however, such methods 
could only deal with idealized continuos potential pro- 
files, and could not include physical limitations in the 
alloy concentrations; for example, the optimized poten- 
tials proposed in Refs. fd,|q cannot be implemented as 
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strictly type-I structures. 

Recently, we have introduced a new numerical strategy 
which allows to optimize the performance of a given de- 
vice with respect to all geometrical parameters and exter- 
nal fields simultaneously^!. In short, our approach uses 
an artificial intelligence tecnique, namely, evolutionary 
programming^, to solve the inverse problem of design- 
ing a material with preset electronic properties, including 
complex non-linear constraints, such as multi-subbands 
transition energies; by exploring the multi-dimensional 
parameter space, the algorithm uses a genetic paradigm 
to find the regions providing the best performances, and 
focuses the search in those regions. Using a population- 
based strategy, and starting with a completely random 
initial set of potential profiles, the algorithm efficiently 
relaxes toward one or more optimal devices without lim- 
iting the search to a given class of structures. Details 
and demonstration of the algorithm can be found in Ref. 
|l9| . The efficiency of our algorithm allows for a quan- 
titative modelling of the structure, including, e.g., fi- 
nite band-offsets, space-dependent effective masses, band 
non-parabolicity, and physical limitations on alloy con- 
centrations. 

In this paper we apply our optimization strategy to 
determine the optimal composition for SHG from in- 
tcrsubband transitions in the conduction band of the 
GaAs/AlGaAs class of materials. We report the com- 
position and potential profile for an optimal device with 
zero electric field. We also study optimal devices with 
a finite F, but we find that, in general, devices with a 
non-zero electric field have an estimated SHG which is 
weaker than devices with F = 0. 



DETAILS OF THE CALCULATION 

We consider an heterostructure composed of a large 
number of Al^Gai-^As layers, each consisting of n mono- 
layers (ML); in order to limit ourselves to type-I het- 
erostructures, the Al molar fraction, x, can take continu- 
ous values in the range [0,0.4]; the GaAs/Al .4Gao.6As 
band offset is large enough to accomodate the three 
bound states which are required to give rise to double- 
resonant SHG. The structural parameters x and n of each 
layer and, in a sub-set of runs, also the eletric field Fc3, 
were subjected to the evolutionary dynamics in order to 
maximize \x = |^i2A*23A*3i| at the frequency of a CO2 laser 
(Tlui = 0.116 eV). As mentioned above, in the double- 
resonant regime the intersubband gaps match the radi- 
ation field; the possibility to implement such non-linear 
contraints is a distinct advantage of our numerical strat- 
egy. Since inhomogenous broadening lowers the peaks 
in the response function, our algorithm allows the struc- 
tures to satisfy the resonance constraint only approxi- 
mately; in the language of Refs. |l9|,|o|, this is obtained 
using a "penalty function" with a Lorentzian shape and 
a width r independent of the transition. 



Calculations were performed within a one-band enve- 
lope function description of the conduction band single- 
particle states in the low-density limit. The single- 
particle Hamiltonian, with a step-like potential deter- 
mined by the sequence of Al concentrations^! appear- 
ing at each step of the simulation, was represented in a 
plane-wave basis set normalized to a periodic box; cap- 
ping layers of typically 40 MLs of Alo.4Gao.6As terminate 
each strucure to ensure that the wavefunctions vanish at 
the boundaries. The space-dependent effective mass was 
linearly interpolated between the. GaAs and AlAs masses 
for arbitrary Al concentrations^. Non-parabolic energy 
dispersion is known to have a sizable effect on subband 
energies and wavefunction localization for large confine- 
ment energies; this is of particular concern in the design 
of a device for SHG, where the three bound states in- 
volved span a large energy range above the conduction 
band edge. By explicit calculations, of the type described 
below, we have found that the optimal SHG intensiy can 
be overestimated by ~ 30% if non-parabolicity is ne- 
glected; in our calculations non-parabolicity effects were 
takea-into account using an energy-dependent effective 
massed. 



NUMERICAL RESULTS 

The potential profiles and subband charge densities of 
an optimal structure, labelled SI, obtained with F = 0, 
and a sub-optimal structure, labelled S2, with F ^ 0, 
are shown in Fig. 1. The corresponding compositional 
parameters are reported in Tab. 1. As can be seen from 
the figure, SI has a remarkably simple shape, with sharp 
variations of the Al concentration rather than a graded 
composition (indeed, sub-optimal structures with a more 
graded composition were obtained during the simula- 
tions; see also S2). Structure SI is remarkably similar 
to the step-graded device considered in Ref. |[ but a) 
with a sharp (high and narrow) barrier between the two 
steps, and b) with a graded region sideways to the deeper 
well, which allows the charge density of the third sub- 
band to slightly delocalize out of the well region; from 
a large set of runs, we conclude that both these features 
are characteristics of the best performing structures. The 
calculated optimal value [i = 3.27 nm 3 of SI exeeds by 
37 % the value /i = 2.39 nm 3 estimated in Ref. ^. 

The structure S2 is a typical example of a sub-optimal 
device with an optimized electric field substantially dif- 
ferent from zero, which appeared in our simulations. The 
structure has a graded double quantum well shape; again, 
a graded region is present sideways to the largest well. As 
shown in Tab. 1, the estimated value /i = 3.126 nm at 
the optimal field F — 43kV/cm, which is ~ 4% smaller 
than for SI. 

Indeed, calculations show that, in general, an electric 
field does not improve SHG if compositional parameters 
are optimized. In Fig. 2 we show the value of \x for 200 
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structures which are optimized with respect to structural 
parameters and the electric field (open dots); the devices 
SI and S2 are also indicated (solid triangles). As it is 
apparent, the best performing structures tend to have a 
small or negligible field; the possibility to achieve large 
SHG is substantially decreased for fields larger than ap- 
proximately 60 kV/cm. 

Although the best performing device, SI, is obtained 
for a null electric field, an optimally designed biased het- 
erostructure,may be desirable if field-controlled SHG is to 
be achieved!!. For example, S2 is designed to be doubly- 
resonant at F = 43kV/cm; when the field is switched 
off, the intersubband gaps become 0.132 eV (1 — > 2) and 
0.103 eV (2 — > 3), i.e., the resonance condition is lost; 
accordingly, the SHG drops by two orders of magnitude, 
although the value of /i is barely affected. 

The benefit from huge non-linearities which can be ob- 
tained by quantum confinement can be limited by ab- 
sorption from the material at the SHG frequencies; it 
has been suggestedQ that optimal conversion efficiency 
might be obtained with /Z12/ V13 — 1- O n the other 
hand, our simulations show that optimal quantum (i.e., 
confinement-induced) SHG efficiency is generally ob- 
tained with ^12/^*13 — 2 (this is, e.g., the case for devices 
SI and S2), which reduces the conversion efficiency by a 
factor of ~ 2. 

In summary, the present investigation of the optimal 
composition for double-resonant SHG in GaAs-based de- 
vices suggests that the largest quantum confinement- 
induced SHG is obtained in asymmetric double-quantum 
wells with narrow and high barriers and a graded com- 
position sideways to the largest well. Furthermore, an 
electric field does not necessarily improve the SHG if 
structural parameters are optimized. 

I acknowledge several useful discussions with Fausto 
Rossi. 
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TABLE I. Compositional parameters for the optimized devices SI and S2. SI is obtained with F = 0, while S2 is obtained 
including F in the optimization. For each structure, the width of each layer, n (in ML) , and the corresponding Al concentration, 
x is reported; two capping layers with x = 0.4 enclose each structure. Inter-subband dipole moments, fiij, and /i — 1^12/^23^31 [ 
are also reported in the table. Calculations have been performed with Tllo — 0.116 eV and T = 7meV; the resulting inter-subband 
gaps match %uj within 1 meV. The correponding potential profiles and subband charge densities are shown in Fig. 1. 
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FIG. 1. Potential profiles of the optimized devices SI and S2; the corresponding compositions are detailed in Tab. 1. Also 
shown are the charge densities of the three bound states involved in the SHG, vertically shifted by the confinement energy. 
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FIG. 2. Dipole moment fi vs electric field for a set of 200 devices (open dots) which have been optimized by evolutionary 
dynamics with respect to compositional parameters and electric field. Also shown are the values corresponding to devices SI 
and S2 (full triangles). 
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